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The existence of a tachocline in the Sun has been proven by helioseismology. It is un- 
known whether a similar shear layer, widely regarded as the seat of magnetic dynamo 
action, also exists in a protoneutron star. Sudden jumps in magnetic diffusivity rj and 
turbulent vorticity a, for example at the interface between the neutron-finger and 
convective zones, are known to be capable of enhancing mean-field dynamo effects 
in a protoneutron star. Here we apply the well-known, plane-parallel, MacGregor- 
Charbonneau analysis of the Solar interfacial dynamo to the protoneutron star prob- 
lem and calculate the growth rate analytically under a range of conditions. It is shown 
that, like the Solar dynamo, it is impossible to achieve self-sustained growth if the dis- 
continuities in a, rj, and shear are coincident and the magnetic diffusivity is isotropic. 
In contrast, when the jumps in 77 and a are situated away from the shear layer, self- 
sustained growth is possible for P < 49.8 ms (if the velocity shear is located at 0.3i?) 
or P < 83.6 ms (if the velocity shear is located at 0.6R). This translates into stronger 
shear and/or a-effect than in the Sun. Self-sustained growth is also possible if the 
magnetic diffusivity if anisotropic, through the f2 x J effect, even when the a, 77, and 
shear discontinuities are coincident. 
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1 INTRODUCTION 



The origin of neutron star magnetic fields is still unclear, especially the impressively high magnetar fields. There are two 
popular scenarios: (1) fossil, where the progenitor star's mag netic field is a mplified by compression during core collapse 
and is frozen into the highly conduct ing compact remnant I Woltier 19641 : Rudermanl 11972 : Braithwaite fc Spruit] 12004; 
Ferrario fc Wic kramasinghc 2006, 2008 ); and ( 2) dynamo, where the seed field is amplified by a convective dynamo in the pro- 



toneu tron star jRuderman fc Sutherlandll973l ; Duncan fc Thompson 1992; T hompson fc Duncanll993l : iBonanno. Urpin. fc Belvedere! 
2005). These origins are not mutually exclusive; for example, remnant fossil fields can survive and be amplified by dynamo 



processes afterwards. 

Neutron star dynamo theory often draws concepts and analogies from Solar dynamo theory. In the context of the Solar 
dynamo, observations from helioseismology show that there is a layer of maximum shear, located between the radiative 
zone (which rotates as a solid body) and the convective zone (which rotates differentially, slower at th e poles than at the 



equator). Called the tachocline, thi s shear layer is centred at 0.69 with a thickness of 0.09 (|Goode et al 



1991 



Tomczvk. Schou. fc Thompson 1995 : Kosovichev 1 19961 : iDikpati fc Charbonneau 19991 : Christensen-Dalsgaard fc Thompson 
2007! ) . It is assumed to be where the toroidal field con ce ntrates: magnetic flux in the convective zone is expelled on time-scales 



of da ys to months ( Parkei 19751 : Glatzmaier 1985al lbl: iFerriz-Mas 19961 : Dikpati fc Charbonneau 19991 : Weiss fc Thompson! 



l2009h . While the situation is less clear in protoneut ron stars, steep angular velocity grad ients are observed in high-resolution 
numerical simulations of core-collapse supernovae (|Ott et al.ll2006l : iBurrows et al.l 120071 ). Moreover, jumps in the magnetic 
diffusivity 77 and turbulent vorticity a between layers, cause d by different species of instabil i ties (convective in the inner region , 
neutron finger in the outer region), have been hypothesized (jMiralles. Pons, fc Urpinll2000l :l IBonanno. Rezzolla. fc UrpinlbooA 
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Table 1. Summary of the dynamo cases discussed in this paper 



Case 



[rj\, [a], [VO] coincident 
[77], [a], [Vfi] not coincident 



Sec. 2.1 Sec. 
Sec. 2.2 



Bonanno. Urpin. fc Belvedere! 20051 ). Sudden jumps in Tj and a enhance mean-field dynamo effects in a protoneutron star: 
poloi dal fields of ~ 10 9 T can be generated in a protoneutron star rotating with period of 0.1 s IjBonanno. Urpin. fc Belvedere! 
20051 ). 

In this paper, we calculate the growth rates of a dynamo operating at a boundary where the magnetic diffusivity and/or 
shear jump discontinuously. We base our analysis on the Solar interfacial dynamo model proposed bv lMacGregor fc Charbonneau 
1 19971 ) but apply it to a protoneutron star, where ma gnetic diffusivity and turbulent vort icity are expected to be discontinu- 
ous between the convective and neutron- finger zones I Bonanno. Urpin. fc Belvedere! 2005 ). Furthermore, the scalar magnetic 
diffusivity 77, which appears on the right-hand side of the magnetic induction equation is actually an anisotropic tensor r\ij 
in general. When analysing the interfacial dynamo, we allow for the possibility of an anisotropic rjij (i.e., one containing off- 
diagonal elements), allowing inte ractions between the homogenous and isotropic turbulence and the inhomogene ous magnetic 
field, known as the fi x J effect ( Rogachevskii fc Kleeorin 20031 ; iRudiger fc Hollerbach 2004 ; Gressel et al. 20081 ) . We wish to 
see if an interfacial dynamo with the above features can operate in a protoneutron star. 



The structure of this paper is as follows. Sec. 2.1 applies the interfacial dynamo model of iMacGregor fc Charbonneau 



1 19971 ). with coincident discontinuities in 77, a, and velocity shear Vf2, to a protoneutron star. Sec. 2.2 extends the analysis to 
the situation where the discontinuities in 77 and a do not coincide with the shear layer. In Sec. 3, we return to the situation 
where the discontinuities in 77, a, and Vfi coincide, but generalise it to include the fl x J effect. The general, non-coincident, 
tt x J case is not treated in this paper, as s umm arized in Table 1. We discuss certain problems with our calculations, and 
indeed those of IMacGregor fc Charbonneau! (|1997r l. in Sec. 4. Lastly, we summarize our findings in Sec. 5. 



2 INTERFACIAL DYNAMO IN A PROTONEUTRON STAR 
2.1 MacGregor &: Charbonneau (1997) interfacial dynamo 

The interfacial dynamo of Parker ( 19931 ) and MacGregor fc Charbonneau ( 1997 ) sits at the tachocline, where shear generates 
toroidal field, which diffuses to the layer where turbulent vorticity (the 'a-effect') converts it into poloidal field, which in turn 
converts bac k into toroidal field via the shear (in this model, both shear and a-effect are sharply localised near the tachocline). 
As shown bv lMacGregor fc Charbonneau! (|1997l ). this dynamo can be stable and self-sustaining in the Sun. We investigate if 
a similar process, operating at the convective-neutron finger transition layer in a protoneutron star, can also be stable. 

The idea of the tachocline comes from Solar physics. The transition between the (inner) radiative zone (where energy 

generated by fusion reactions in the core travels outwards by radiative transport) and the (out er) convection zone (where energy 

moves by convect i on) occurs at 0.71 Rp, where Rp is the Solar radius, independent of latitude (jChristensen-Dalsgaard, Gough. fc Thompson 



199ll ; lBasulll99g| ; iBasu fc Anti3l200lt [Jones. Thompson, fc Tobiasl hoosl ; IPetrovav fc Christensenlboidl ). The radiative zone 



rotates as a solid body, while the convective zone rotates slower at the poles than the equator. The layer of maximum 
shear flow between the radiative and the convective zones, which is the tachocline, is centre d at 0.69 Rp } with a thick- 



ness of 0.09 Rp dGoode et all 1991 


; Tomczvk. Schou. & Thompson 1995 


Christensen-Dalsgaard & Thompson 


2007). 



2005) 



The mean- field induction equation in a turbulent protoneutron star can be written as ( Bonanno. Urpin. fc Belvedere! 



dB 

dt 



= V x (v x B + aB) - V x (77V x B) 



(1) 



where B is the magnetic field, v is the mean fluid velocity, 77 is the turbulent mag netic diffusivit y, and a, the so-called 
'a-parameter', is a pseudo-scalar measuring the efficiency of the dynamo. According to iParker jl993l ). dynamo processes are 
concentrated around the tachocline, where a is sourced by local conv ective motion, which can vary across the tachocline 
( Parkenfliiil ; Izhang. Liao. fc Schubert] |2004| ; I Weiss fc Thompson! [2009! ) . The introduction of discontinuities in a, v, and 77 
supports solutions to Eq. (1) in the form of exponentially growing travell ing waves at the interface, creating a source term 
for the poloidal component and the possibility of a self-sustaining dynamo ( Parker 19931 : Dikpati fc GilmanlEoOll ). This is the 
basic outline of the so-called a-il dynamo mechanism: t he convective fluid m otions give rise to a source term for the poloidal 
field and the a-effect generates the toroidal component 1 Charbonneau 20051 ) . 

If the velocity shear (which generates the toroidal field component) and the turbulent a-effect (which generates the 
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poloidal field component) are co -located, the convective fluid motions tend to suppress the a-effect and resistive diffusion, 
thereby dampening the dynamo 1 Parker 19931 ; MacGregor fc Charbonneau 1997). If the a-Q. model is modified by separating 
the shear layer and the a-effect by some distance, the suppression is lifted and the ratio of the azimuthal fields (which act 



as sources o f the poloida l fields) in the top and bottom layers varies inversely as the square r oot of the ratio of magnetic 
diffusivities (jParkerll 19931 : iMacGreeor fc CharbonneaJ 1997h . iMacGreeor fe Charbonneau! (1 1997 ). 



motiv ated by observations 

regarding the thinness of the tachocline |Kosovkhev[1l99rl ). modified the aforementioned Parker l| 19931 ) dynamo further by 
introducing a sharp discontinuity in the velocity shear (as well as a, and magnetic diffusivity rj). For this interfacial dynamo 
model, where the dyna mo operates in the region around the discontinuities, they found growing modes confined to finite 
ranges of wavenumber k IjMacGregor fc Charbonneaulll997V ). 

To model the thin region around the tachocline, we adopt a semi-infinite Cartesian geometry: the poloidal direction is 
taken to be in the x- and z-directions, the toroidal direction is in the y-direction. The field B is assumed to be symmetric in 
y, viz., 



B = B x (x, z,t)e x + B(x,z,t)e y + B z (x, z,t)e z ; 
Its poloidal components can be derived from one quantity, the vector potential A = A(x, z,t)e y : 



(2) 



dA dA 

B x e x + B z e z = V x A(x,z,t) = - — e x + — e z 

dz ox 



(3) 



Following iMacGregor fe Charbonneau! (jl997j), we take the large-scale velocity field v to be constant and pointing in the 
y-direction, with a nonzero vertical gradient, viz. v = v(z)e y . Equation (1) now splits into poloidal and toroidal components: 



dA 
~dt 



= V 



it 

dx 



d 2 

2 + o^) A + aB > 



(4) 



dB 
~dt 



' dx 2 dz 2 



B+-— + dr) dB 

dz dx dz dz 



(5) 



Equation (4) follows from both the x- and z- components of Eq. (1) , which are identical after (but no t before) integration 
with respect to z and x respectively (this issue is revisited in Sec. 4). MacGregor fe Charbonneau! 1 1997 ) chose 



rj = rjx + (r]2 - -ni)Q(z - h), 



(6) 



a = 7o<5(z - d), 



dv 
dz 



= LUo5(z), 



(7) 
(8) 



where the subscript '1' ('2') indicates the lower (upper) region, and O(z) and S(z) are the Heaviside and Dirac delta functions 
respectively. In other words, the discontinuities in a and rj occur at z = d and z = h respectively, while the discontinuity in 
velocity shear is located at z = 0. 

When applying this model to a protoneutron star, r/i — rj c is the magnetic diffusivity in the (lower) convective region, 
and rj2 = Tj n f is the magnetic diffusivity in the (upper) neutron-finger unstable region. Physically, the different r\ is because, 
early in a protoneutron star's life (< 100 ms after core collapse), the inner region of the star is conv ectively unstable (driven 
by entropy gradient), whereas the o uter region is neutron- finger unstable (driven by lepton gradient) i Bruenn fc Dineva 199rj ; 
Bonanno. Urpin. fc Belvedere! 2005 ). On the other hand, our of the 5-function for a, to locate the a-effect sharply at the 
interface, is more of a mathematical convenience; we assume that a, which comes from local convective turbulent, changes 
sharply at the interface due to some unspecified fundam ental physical difference between the convective and neutron-finger 
region s I Parker 19931 : MacGregor fc Charbonneau 1997). To make contact with the work of Bonanno. Urpin. fc Belvedere! 
(2005). we consider h — d only. We take h — d — in Sec. 12.21 then investigate h = d / in Sec. 12.31 



2.2 Coincident discontinuities in a, 77, and velocity shear 



Integ rating Eqs. (4) and (5) across the interface, implies the following jump conditions on the field variables ( M acGregor fc Charbonneau 
1993): 



[A] = 0, 

[B] = 0, 



(9) 
(10) 
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Vc + Vnf 



'OA' 

dz 



= -70-B, 



as" 

' dz 



= — wo- 



dA (z = 0) 
dx 



(11) 



(12) 



Square brackets denote the differences between the values taken at z = +e and 2 = — e (where e is a small number). These 
conditions are obtaining by requiring that the normal and tangential fields are continuous across the interface [mathematically, 
this ensures that there are no derivatives of infinite quantities in Eqs. (4) and (5)] and by integrating Eqs. (4) and (5) across 
the interface. We discuss the derivation of these boundary conditions in more detail in Sec. 4. 

The solutions to the linear equations (4)-(5) away from the interface (i.e., z 7^ 0, a — 0, dv/dz = 0, and 77 = 771,2) can 
be written as A, B oc exp(±q Cin /z) exp[crf + i(kx + OJt)], where q% n t = fc 2 + (a + iw)/r] c>n f, after discarding the exponentially 
growing solutions at z — > ±00. Upon substituti ng the exponential trial solutions f or A and B into the jump conditions Eqs. 
(9)-(12), the dispersion relation takes the form IjMacGregor fc Charbonneaulll997l ) 



l + n(l + 2s) + (n + 1)[(1 + s)(l + ns] 
with n = Vc/Vnf, s = (a/r] c k 2 ) + i(uj/rj c k 2 ), and 



11/2 



■iv = 0, 



(13) 



V = 2^070 /[Vnf(Vc + 7?n/)fc]- (14) 

We now calculate the growth rate and frequency of the d ynamo wave solution given by Eq. (11 3 ^ for a protoneutron star. 
Typ ical magnetic diffusivities are r\ a g 10?7 n / ~ 10 8 m 2 s _1 i Bonanno. Urpin. fc Belvedere! 20051 '). so we take n — 10; note 
that lMacGregor fc Charbonneaul (1997) only considered n < 1. In Fig. [T] we plot o-/?7 c and oj/t/'c versus k (in units of m~ ) for 
v — 20 (solid curve), !/ = 10 (dashed curve), v — 5 (dashed-dotted curve), and = 1 (dotted cur ve). We find that the dynamo 



never grows (a is always negative), rega rdless of v or fc. This is consistent with earlier findings i Deinzer fc Stix 197ll : Parker 



1993; iMacGreeor fc Charbonneaul 1 19971 ) that, when the a-effect and velocity shear are at the same location, the convective 
fluid motions tend to suppress the a-effect and the dynamo amplification: th e a-effect does not have enough time to generat e 
the toroidal component of the field before being overcome by velocity shear ( Parker 19931 ; MacGregor fc Charbonneau 1997 ). 

One might wonder whether a sufficiently large v could trigger a growing mode. This is never the case. To see this, note 
that the real part of s which satisfies the dispersion relation Eq. (|13|) can be written as 



Re(s) = 



-{n + 1) ± \/(n- l) 2 -Anv 2 /(n- l) 2 



2n 



since n and v are always real and positive. Growth requires 

(n-1) 



> (n + 1) , 



(n-1) 2 



(15) 



(16) 



which is never satisfied for any real value of v. 



2.3 Coincident discontinuities in a, 77, away from the discontinuity in velocity shear 



We now separate the shear layer and a discontinuity, as lMacGregor fc Charbonneaul (|l997T ) did for the Solar dynamo in an 
effort to circumvent Parker's (1993) quenching mechanism. There are now three zones to investigate: z < 0, < z < d, and 
d < z. For z < and z > d, one physically meaningful mode exists for each of A and B, after discarding the divergent 
solutions oc exp(— q c z) and exp(q n fz) in the regions z < and z > d respectively. For < z < d, two independent modes exist 
for each of A and B. The jump conditions take the form 



rjc + r]„f 



[A] = 0, 

[B] = 0, 
OA 



dz 

dB 

dz 



-70S, 



= 0, 



(17) 
(18) 

(19) 
(20) 



at z = h = d, and 
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Figure 1. Plots of a/r/c (top) and ui/rj c (bottom) versus wavenumber k (where a is the dynamo growth rate, lo is the dynamo frequency, 
and r/c is the magnetic diffusivity in the convective regio n), with the magnetic diffusivity ra tio n = r] c /r) n f = 10 (where r) n f is the 
magnetic diffusivity in the neutron-finger unstable region) IIBonanno. Urpin, k. Belve dere 2005), for the case where there is no separation 
between the discontinuous shear layer and the discontinuities of r\ and o, for four different values of v = 2wo7o/['?n/ (ilc + T]nf)k] (where 
ojo is the discontinuous velocity shear and 70 is the discontinuity in a): v = 20 (solid curve), v = 10 (dashed curve), v = 5 (dashed-dotted 
curve), and u = 1 (dotted curve). 



[A] = 0, 

[B] = 0, 



OA 

dz 



= 0, 



dB 



dz 



dA(z = 0) 
dx ' 



(21) 
(22) 

(23) 
(24) 



at z = 0. We therefore have eight equations t o solve for eight unknowns. 

We expect to see growth in this situation: iMacGregor fc Charbonneaul |l997 ) found a band of dimensionless wavenumbers 
K — 2kd which admits growing modes, as seen in their Fig. 2. The width and centre of the band depend on n and v through 
the dispersion relation 



1 + n{l + 2s) + (n + 1)[(1 + s)(l + ns)] 



1/2 



iv exp[ 



-n(l + s) 1/2 } =0. 



(25) 



Note that it is possible for the velocity shear discontinuity to be located in the upper region (the neutron finger region) 
instead. To derive the analogue to Eq.[25]for this particular case, one simply replaces r} c with rjnf in Eq. (24), h with —h, and 
d with —d everywhere. We only consider the case where the velocity shear discontinuity is located in the lower region (the 
convective region) in this paper, because, as argued by Bonanno, Rezzolla, & Urpin (2003) and Bonanno, Urpin, & Belvedere 
(2005), this region is expected to be more convectively unstable. 
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Taking n = 10 again, we plot the real (representing growing/decaying modes, the four lower black curves) and imaginary 
(indicating the presence of dynamo waves, the four top blue curves) parts of s in Fig. [2] as functions of k for v = 50 (solid 
curves), v — 30 (dashed curves), v — 20 (dashed-dotted curves), and v — 10 (dotted curves). We find that the real part of s 
decreases for k > 1. Why is this so? The toroidal field B (generated by shear at z = 0) diffuses to the upper interface at z = d, 
wher e it is converted into the poloidal c omponent A by the a-effect, and vice versa, creating a closed, self-sustaining dynamo 
loop ijMacGregor fc Charbonneaulll997h . For the dynamo to be self-sustaining, one needs source terms for both the toroidal 
component (B) and the poloidal component (A) in Eqs. (4)-(5): both A and B need to grow and sustain each other. As k 
increases, the first terms on the right hand side of Eqs. (4) and (5), which are sink terms oc — (r/k 2 ), become more negative 
and cannot be balanced by the source terms a or du/dz. When the dissipation time-scale (rjk 2 ) -1 is shorter than the diffusion 
time between layers (~ d 2 /rj), waves with higher k damp faster than they can diffuse and the dynamo processes in Eqs. (4) 
and (5) cannot regenerate the field components any longer, regardless of the shear or the strength of the a-effect. In other 
words, the shear layer and the a-effect layer cannot be co-located (as shown in Sec. 2.2), but they cannot be too far apart 
either. 

We also see that none of the tested values of v show positive Re(s) within the n domain in Fig. [2l i.e. none o f the c ases 
exhibit a growing dynamo within 0.01 ^ k ^ 100. This is in contrast to the results of MacGreeor fc Charbonneaul (1997); as 
we see in their Fig. 2, v ^ 20 gives a growing dynamo at k — 1. However, if one increases v sufficiently, one does eventually 
achieve Re(s) 0. To see the relation between the physical parameters uo, 70, r\, and k more clearly, let us define 



Q = 4uJo-yod/rj n f(rj c +??„/), (26) 
i.e., Q — vk, and let Q c denote the minimum value of Q (at some fixed n) required to achieve a growing dynamo. In the 



specific case n = 10, expected in a protoneutron star (jBonanno. Urpin. fc Belvedere! 120051 ). we find empirically 



with c = 4.485 x 10 2 , ci = 0.553, c 2 



Q c w c k C1 exp(c2« ' +C3K 1 '), (27) 

-.391, and C3 = —9.187 to within ~ 10 per cent accuracy. In Fig. [3] we plot Q c as a 
function of k, both the true values (dots) an d the approximation (solid curve). We gener ally find Q c values larger than thos e 
found bv lMacGregor fc Charbonneaul (|1997l ). This is because we have, in the notation of lMacGregor fc Charbonneaul (|1997l ). 
771 > 772, i.e. the waves dampen faster (than in the Solar case) in the layer between the shear interface and the a interface, 
so the source terms (oc 70 an d ujq) must be larger to compensate . In other words, in the protoneutron star regime, where 
rj c is expected to be ~ 10r;„/ ( Bonanno. Urpin. fc Belvedere 2005), larger amounts of shear (alp ) and/or turbulent vorticity 
discon tinuity (70) are needed for an operational dynamo, compared to the Solar regime studied by MacGreeor fc Charbonneaul 
|l997h . 

Let us examine this protoneutron star example further. Taking r\ c = 10?7„/ = 10 s m 2 s _1 and R — 10 4 m as the radius 
of the protoneutron star, we find that the dynamo grows for 



10 5 m 



7(1 



10 7 m 2 s" 1 



10 3 



10 7 m 2 s" 1 



C K min eX P( C 2« 



1/2 



_l/3 . 



«: 3.636, 



(28) 



with K m in = 2fe m i n d, where fc m i n is the smallest wavenumber. This wavenumber is & m i n = R^ 1 , where R c is the distance of the 
shear layer from the centre of the protoneutron star. We can then estimate ljo as the tangential velocity of the rotating star 
at R c , i.e. ujo « 2nR c /P (this is, of course, the upper limit of the velocity shear). If the thickn ess of the discontinuity layer is 
A.B = 0.025.R and a w ML « 6 x 10 3 tt/P, where L w 3 x 10 3 m is the pressure length-scale ijBonanno. Urpin. fc Belvedere! 
2005), then 70 ~ aAR follows from Eq. (7) (i.e., the 'discontinuous' jump in a takes place over AR). Substituting these 
parameters into the inequality (|28[) and taking d = AR, we find that the dynamo grows for P ^ 49.8 ms, for R c = 0.3ii, 
or P ^ 83.6 ms, for R c = 0.6i?. Additionally, one can now calculate the coefficients q c ,nf ( as defined in Sec. 2.2), which 
determine the radial length-scale of the dynamo-wave magnetic field, for the above protoneutron star example. When the 
dynamo is critical [i.e., when Re(s)=0, when inequality (28) is satisfied] for k find \q 2 \ = 3.24 x 10 -7 m" 2 and 

\ql f \ = 3.08 x 10~ 6 m -2 for R c = 0.3R and \q 2 \ = 1.24 x 10~ 7 m -2 and \q 2 f \ = 1.21 x 10~ 6 m~ 2 for R c = 0.6R. Hence, the 
dynamo fields penetrate many multiples of the transition layer thickness into the convective and neutron finger zones. 

We find that the dynamo is more easily excited when the interface is farther away from the stellar core. The reason for 
this is obvious: the farther the interface ; the greater the shear, because we ta ke the shear simply to be equal to the tangential 
velocity at the interface. In comparison. iBonanno. Urpin. fc Belvedere! (2005) found that their dynamo develops when P ^ 1.4 
s (for R c = 0.3i?) or P ^ 1 s (for R c = 0.6R). Thus, in addition to the higher threshold periods. (Bonanno. Urpin. fc Belvedere! 
! 2005h reported the opposite trend to ours: the dynamo is more feasible when the interface is closer to the core. The reasons 



Unlike the Sun, we currently do not have a method of locating the shear layer or the a-layer observationally in a protoneutron star. 
The values used in this section are for illustrative purpose only. 
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Figure 2. Real (the four lower curves, in black) and imaginary (the four top curves, in blue) parts of s = (a/rj c k 2 ) + i(ui / 'r) c k 2 ) (where a 
is the dynamo growth rate, ui is the dynamo wave frequency, r) c is the magnetic diffusivity in the convective region, and k is the dynamo 
wavenumber) , as functions of k = 2kd (where d is the separation between the shear layer and the rj di scontinuity), with magnetic diffusivit y 
ratio n = 7\c/Vnf = 10 (where r] n f is the magnetic diffusivity in the neutron-finger unstable region) jBonanno, Urpin, fc B elvedere 2005), 
for different values of v = 2aJo7o/[»?n/ (Vc + Vnf)k] (where u>q is the discontinuous velocity shear and 70 is the a discontinuity): v = 50 
(solid curves), v = 30 (dashed curves), v = 20 (dashed-dotted curves), and v = 10 (dotted curves). 



for this discrepancy are unclear, bu t we remind the readers that our approach is fundamentally different to that taken by 
Bonanno. Urpin. fc Belvedere! |2005h ; the latter authors considered the enti re star, while we only look at the reg ion near the 
interface. It is possible that there are effects and back-reactions present in the lBonanno. Urpin. fc Belvedere! 1 20051 ) model that 
are absent in ours. While our calculation is less rigorous, we present it to show how an analytic interfacial Solar dynamo-like 
model can be usefully applied to a protoneutron star with some additions and modifications. As it is, the threshold period 
values we obtain are still larger than the birth periods > 10 ms of typi cal radio pulsars and > 5 ms of the known magnetars 
(inferred from shock velocity and associated supernova remnant radii) ( Vink fc KuiperlEooi ). 

To summarize, in this section, we show that an interfacial dynamo model, where the a-effect is situated at or near a 
discontinuity in magnetic diffusivity and velocity shear, can operate in a protoneutron star spinning sufficiently fast. The 
velocity discontinuity separates the convective (lower) region from the neutron finger (upper) region in our case, cf. the 
radiative (lower) region from the convective (upper) region in the Solar dynamo case (MacGregor & Charbonneau 1997). In 
the protoneutron star, n, the ratio between the magnetic diffusivities in the lower region and the upper region, equals 10 
(Bonanno, Urpin, & Belvedere 2005) [on the other hand, MacGregor & Charbonneau (1997) investigated n < 1]. Despite this, 
we find [as Deinzer & Stix (1971) and MacGregor & Charbonneau (1997) did] that, when the a-effect, the discontinuity in 
r/, and velocity shear are coincident, no self-sustaining dynamo action is possible. We need d 7^ so that a sufficient toroidal 
component can be generated by the shear before being converted back into poloidal by the a-effect. If the a-effect and the 
discontinuity in r/ are not coincident with the velocity shear, a self-sustaining dynamo is possible for certain parameter ranges. 
However, if d gets too large, the dynamo waves damp before they can diffuse from one interface to anothei0. We derive the 
formula for the threshold Q c [Eq. (27)] and plot Q c as a function of k in Fig. [3] We note that, for a given k, our system requires 
higher v (i.e., stronger a-effect and/or greater velocity shear), compared to the Solar case of MacGregor & Charbonneau (1997), 
to establish a self-sustaining dynamo [compare our Fig. |3]to Figs. 2 and 4 of MacGregor & Charb onneau (1997)1- For example, 
we see in Fig. [3]that, for k = 1, we need v > 195 for a growing dynamo, whereas, from Fig. 2 of lMacGregor fc Charbonneau 
( 19971 ). we see that the Solar dynamo needs v > 20. 

One should note that, throughout this section (and, indeed, throughout the paper), we model the transition layer between 
the neutron finger a nd convective regions as being infi nitesimally thin for simplicity, following the treatment of the solar 
dynamo interfaces bv lMacGregor fc Charbonneau! (|l997l ). This transition layer can, and perhaps should, be treated as having 



2 For the protoneutron star case presented here, we find that the dynamo becomes ineffective only at d ~ R 
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1000 




Figure 3. The critical Q = vk versus k = 2kd (where k is the wavenumber and d is the separation between the shear layer and the r\ 
discontinuity). Q c is the smallest Q for which there is a growing mode [Rc(s) > 0] , with the ratio of magnetic difTusivities n = rj c /r] n f = 10. 



a small but nonzero thickness [ ~ 0.025/?, according; to Bonanno. Urpin. fe Belvedere! ( 2005 )]. Such a calculation has been 
performed for the solar dynamo: Izhang. Liao. fc Schubert ( 2004 ) treated the tachocline layer as having a nonzero thickness 
[so the velocity gradient dv/dz, which we take to be a discontinuous jump in Eq. (8), is a continuous change over a thin 
region between the radiative and convective zones]. They found that, for cases where the thickness of the tachocline is less 
than 0.1-R, the results are la rgely insensitive to the actual thickn ess chosen [recall that the protoneutron star transition layer 
is expected to be ~ 0.025i? 1 Bonanno. Urpin. fc Belvedere! 20051 )]. It would be worthwhile to verify if thi s is st ill indeed the 
case in protoneutron star dynamos, but the calculation is rather cumbersome Zhang. Liao. fc Schubertl (|200J) had to deal 
with 11 linear simultaneous equations] and lies outside the scope of this paper. As we see above, |q c ,n/| -1 are greater than 
the thickness of the discontinuity layer, thus the layer can be assumed to be infinitesimally thin to a good approximation. 



3 (IxJ EFFECT IN AN INTERFACIAL DYNAMO 

In this section, we explore the possibility of allowing anisotropic jjy (i.e., one containing off-diagonal elements) in the induction 
equation, allowing interacti ons between the homogenous and isotropic turbulence and the inhomogeneous magnetic field, 
known as the flxJ effect (Rogachev skii fc KleeorinlEool : iRiidieer fc Hollerbachl Eool : iGressel et al.l 120081 ) . Turbulence in 
rotating fluids often exhibits a preferred di r ection, set by the rotation axis, independent of other factors like density or 
magnetic diffusivity gradients ( Radler 19691 : Roberts! 19721 : Rogachevskii fc Kleeorin 20031 ). The contribution to the time- 
averaged electromotive force (v x B) from this anisotropic eff ect is linear in gradients of B to leading order , hence the term 
'J7 x J effect', where ft is the local vorticity of the mean flow (|Roberts!ll973 : lRogachevskii fc Kleeorinlboost ). The additional 
electromotive force is generally absorbed into the m agnetic diffusivity tensor rjij in Eq. (1) as off-diagonal, symmetric entries 



I Rogachevskii fc Kleeorin! 



2003 



Gressel et al 



2008). In this section, we discuss the effects of including off-diagonal rjij terms 



in the interfacial lMa(^GregoT*fc^CTiai^oTn^ai' (|l997h dynamo, as applied to a protoneutron star. We ask the question: does 
the inclusion of off-diagonal rjij entries enhance or suppress the growth of the interfacial dynamo, discussed in Sec. 2? 



Firstly, we augment the induction equation by adding an extra term 



c^B 

— = Vx(vxB + aB) 



V x (77V x B) + V x E, 



(29) 



with 



E, 



-VijtjkiVkBi, 



(30) 
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where rj is the isotropic magnetic diffusivity (i.e., the trace of rjij) 
only, which we assume, for simplicity, to be given by 



and rj ij is a tensor containing the off-diagonal elements 




(31) 



The form (29) with (30) and (31 ) is standard in stellar magnetohydrodynamics 1 Rogachevskii fc Kleeorinl2003l ; Riidiger fc Kitchatinov 
20061 ; IBrandenburg et alll2008h . 

Next, we investigate how the new terms modify the interfacial dynamo system in Cartesian geometry. We write the 
magnetic field B as in Eqs. (2)-(3). The poloidal and toroidal components evolve according to Eqs. (4)— (5)] 



dA 
~dt 



= V 



it 

dx 



- ^ dz 2 



dB 
~dt 



= V 



dx 2 



A + olB + fl 



d 



dB 



dvdA dr) dB 
dz dx dz dz dz 



: 2 ^ dz 2 



A 



(32) 



(33) 



with 77, a, and dv/dz defined in Eqs. (6)-(8). Importantly, Eqs. (32) and (33) are obtained directly from MacGregor & 
Charbonneau's (1997) equat i ons (6 ) and (7) with the ji terms added. We argue below, in Sec. 4, that the approach taken by 
MacGregor fc Charbonneaul 1 19971 ) (and, hence, also here) is mathematically improper, but we take it at face value for now 

and proceed by handling the fi terms in the same way. 

What form does (i take? An analytical time- and space-dependent expression h as been derived by Urpir] d2002 ) for the 



a-Q. dynamo. In some cases, it can be calculated nu merically for specific sy s tems ( Brandenburg! 2005 ; Brandenburg et al 



2008; Rheinhardt fc Brandenburg! 201(3 ). According to Riidiger fc Kitchatinov 1 2006h . jx must be nonuniform in the direction 
normal to the plane of the shear (i.e., dfi/dz ^ in our formulation). Hence, for reasons which will become readily apparent 
later when we discuss mathematical subtleties in Sec. 4, we assume here a very simple form for fi, which is time-independent 
and sharply localised, i.e. 



ft = XoS(z). (34) 

Note that [i now has the same form as a, as defined in Eq. (7). This assumption is physically defensible; both a and fi 
ultimately stem from the turbulent vorticity, so if the isotropic a-effect is concentrated in a thin layer (by the 5-function), its 
anisotropic counterpart Tj.j is probably similarly concentrated and takes the same form. First-principles numerical calculations 
of the coefficients \o in a protoneutron star lie outside the scope of this paper. 

In this section, for clarity, we consider only the case of h = d — 0, i.e. where the shear layer coincides with the a and rj 
discontinuities. Without the fi x J effect, this situation is shown to be unstable, with purely decaying modes (see Sec. l2.2[) . We 
ask if the fi x J effect drives growing modes in the system. The boundary conditions, obtained by integrating Eqs. (32)-(33) 
with respect to z across the boundary [and by requiring that no terms in Eqs. (32)-(33) are derivatives of an infinite quantity], 
reduce to [cf. Eqs. (9)-(12)] 



dB 
? 97 



Vc + Vnf 



[A] = 0, 

[B] = 0, 
dA 



dz 



= -70-B, 



= — wo- 



dA(z = 0) d 2 A(z = 0) 



d.i 



(35) 
(36) 

(37) 
(38) 



Eq. (37) is obtained by integrating Eq. (32) with respect to z. It is identical to Eq. (11) (but see Sec. 4). Eq. (38), the 
counterpart of Eq. (12), is obtained by integrating Eq. (33) with respect to z. The analysis of Sec. 12.21 remains valid away 
from the interface. We find A, B oc exp(±q Cjn /z) exp[<rt + i(kx + wt)] as before, with <i£ n / = k 2 + (a + iu))/r) c>n f. Substitution 
of the trial solutions for A and B into the boundary conditions yields the following dispersion relation [cf. Eq. (13)]: 



with 



1 +n(l + 2s) + (n + 1)[(1 + s)(l + ns)] 1/2 - iv + v x = 0, 



(39) 
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— , (40) 

Vnf 



(41) 
(42) 



rick 2 \r] c k 2 
2^o7o 



ui = 2 7o xo 3 

VnfWc + TInf) 

As above, subscripts 'c' Qnf) denote the lower, convective (upper, neutron ringer) region. The parameters n, u, and v\ are 
all real. 

The extra term in Eq. (39), arising from fi, is real and positive, suggesting that solutions with Re(s) > may now be 
possible. By directly evaluating the roots of Eq. (39), we search for the minimum v\ requi red to obtain growing modes, denote d 
as ^i, m in, versus v, for a typical protoneutron star with magnetic diffusivity ratio n = 10 |Bonanno. Urpin. fc Belvederejliooi - ). 



The result is shown as the solid curve in Fig. [3] i^i, m in first grows quickly (by about ~ 100 times) for v < 20, then levels off 
to approach a constant value z/i, m i n — > 8.1 for v > 20. Thus, we find that there is always some v\ for which Re(s) > 0; in 
other words, one can always find some v\ (or, equivalently, xo) which leads to growing modes for all n. For completeness, we 
also plot fi,min for a smaller magnetic diffusivity ratio n = 3, as the dashed curve in Fig. [4] We find that V\ >m i n approaches a 
smaller value fi, m i n — > 0.62, for v > 10, but the conclusion stands. 

We designate the value approached by vi, m m as v — > oo as vi t t (i.e., vi t t is the minimum value of v\ needed to guarantee 
growing modes for a given n, regardless of v) and we plot vi t t versus n in Fig. [5] We find that v\x peaks at n ~ 250. Note 
that this behaviour is essentially independent of v or any other physical parameters of the system. At the time of writing, the 
physical justification of this behaviour remains unclear. 

By way of illustration, we evaluate the condition for growth for a typical protoneutron star, with r\ c — 10?7„/, fi,t = 8.131, 
and hence 



( 2* )( *2 )( Hal r 2 > 4 472 (44) 

V10 8 m 2 s- 1 / V10 7 m 2 s- 1 / V10 7 m 2 s- 1 / " ' ' K ' 



Let us assume that the protoneutron star rotates with period 100 ms, so that 70 ~ aAR ~ 1.5-7T x 10 7 m 2 s _1 (Sec. 2.3). 
Then, according to Eq. (43), v\ 8.131 implies xo ^ 9.49 x 10 7 m 2 s" 1 . If we also assume that the jump in fj, occurs in a layer 
of thickness ~ AR, then, by the definition of (1 according to Eq. (34), we can approximate /i ~ xo/AR ~ 3.796 x 10 5 m s _1 . 
This value of \i corresponds to 2.01 a. In summary, for n = 10 and /1 > 2.01a, the f2 x J effect drives a growing dynamo in 
a protoneutron star where the velocity shear a nd a, r\ discontinuities are coincident. In this sit uation, one expects no growth 



without the ft x J effect, as shown in Sec. 12.21 (Parker 19931: iMacGreeor fc Charbonneaulll997l ) 



4 MATHEMATICAL SUBTLETIES 



Sections 2 and 3 of this paper follow closely the analysis of iMacGregor fc Charbonneaul (| 19971 ) . In the process of applying 
and extending the original paper, we encountered some physical and mathematical inconsistencies, which we set out in this 
section. Resolving these matters, which also affect our calculation, is a substantial task lying outside the scope of this paper. 

As a first step, we write out fully all three Cartesian components of the induction equation, Eq. (29), for the field given 
by Eqs. (2)-(3), with a, v, and 77, and fi varying only with z. The result is 



dB 
~dt 



dvdA 
dz dx 



d A 
dzdt 

d_ 

dz 



d(aB) d_ 

dz dz 



— + — ) I 

cte 2 dz 2 



d_ ( dB_ 

dz v dz 



dA 



+ a 



d 2 A 

dx 2 



dB 



d 2 B d_ 

^ dx 2 dz \ dz 



h — 

dz 



/' 



cte 2 



dz 2 



A 



d 2 A 



dB 



d 



dxdt 



dx 



'dx 



d 2 



A 



d 2 B 



dxdz 



(45) 
(46) 
(47) 



v dx 2 dz 2 

Note that the placement of a and r\ inside or outside particular brackets is deliberate and crucial for maintaining full generality 
We now observe the following points. 



The second term on the right-hand side of Eq. (46) is entirely missing from Eq. (7) in the original paper bv lMacGregor fc Charbonneau 



(1997|). This term, arising from the V x (qB) term in the induction equation, is important in view of the discontinuity in a. 
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Figure 4. Minimum fl X J discontinuity Ui jIn ; n versus jump amplitude for the cases n = rj c /r] n f 
curve), where !^i, m j n is the minimum v\ required to obtain growing modes. 
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Figure 5. The minimum Si X J discontinuity required to guarantee growing modes, vi,t, as a function of n = »? c /»7n/ (with jump 
amplitude v — > oo in all cases). 



MacGregor fc Charbonneaul (|1997l ) regarded Eqs. (45) and (47) to be equivalent and hence deemed (45) to be redundant. 



While this assumption is valid away from any discontinuities, complications arise when one tries to derive jump conditions 
across z = 0. Eq. (47) can be safely integrated with respect to x, since everything is continuous in x. In contrast, the derivatives 
with respect to z in Eq. (45) cannot be so easily undone. Equations (45) and (47) differ by terms which jump across the 
interface, like [otB] and [fi(dB/dz)]. In turn, this creates a mathematical dilemma: if there are only two fields A and B, 
how can they be governed by three independent equations Eqs. (45)-(47)? One possible resolution is to require J d;z(45) and 
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J dx(47) to be identically equal, which implies [aB] = and [(j,(dB/dz)] = 0. Because a and \x are Dirac S functions, this in 
turn requires [B] — [dB/dz] = 0, and the problem reduces to the trivial A = B = solution. 

• When attempting the case fx ^ 0, we find that the simplest assumption, that of a uniform and constant fx, yields 
solutions that are linear superpositions of four independent modes in each layer, multiplied by eight unknown coefficients. It is 
unphysical to have any infinite terms in the equation, except for the highest derivatives. For example, a term like d 3 A/dz 3 can 
be singular (e.g., a S function) at z = if d 2 A/dz 2 is discontinuous ([d 2 A/dz 2 ] = 0), because integrating over this singularity 
gives something finite; but it is unacceptable for dA/dz to be discontinuous too. Thus, a term like d 3 A/dz 3 [e.g., the last 
term in Eq. (46) with fi uniform] implies that the dA/dz term in Eq. (46) must be continuous. However, we find that we only 
have seven conditions to solve for eight unknowns (namely, the aforementioned coefficients of the trial solutions for A and B 
in the upper and bottom regions), implying that the problem is badly posed. 

• Suppose instead we try fx = xoS(z), which is justified physically in Sec. 3. We now have only two modes in each layer, 
like for y, = (as opposed to four). However, it is mathematically improper to have terms which are products of infinities. For 
example, a term like y d 2 A/dz 2 is effectively a product of t wo 8 functions, if ji = xoS(z) and dA/dz is discontinuous at z = 0, 



as it is in Eq. (14b) of iMacGregor fc Charbonneaul |l997]). We then find that some of the boundary conditions can only be 
reconciled with each other if r/ is continu ous across the interface, whi c h con tradicts our basic postulate. In Sec. 3, we have 
thus ignored Eq. (45) completely, copying IMacGregor fc Charbonneaul ((1997), as well as the potentially pathological product 
terms, namely the last term of Eq. (46) and the last term of Eq. (47). 

• One can instead let /_t vary continuously between the bottom and top regions, e.g. as a linear ramp. By this expedient, 
one avoids the aforementioned mathematical issues. However, one now has three regions to consider (the top and bottom 
ones plus one transitional region between them) and, therefore, twelve unknown c oefficients. An examp le of this approach at 
work is given by Chandrasekhar in his chapter on the Rayleigh- Taylor instability (|Chandrasekharll 19811 ) . We do not take this 
approach in this paper because our goal is to apply the Parker-MacGregor-Charbonneau dynamo model to a protoneutron 
star and perform direct comparisons. 

The above issues must be resolved before one can obtain self-consistent solutions to the interfacial dynamo problem in 
protoneutron stars and even in Solar contexts. While a resolution is out side the scope of this paper, we hop e that the results 
of Sections 2-3 are still physically indicative, at least to the extent that IMacGregor fc Charbonneaul (1 19971 ) is a fair model of 
the Solar dynamo, and c an be used to guide further calcula tions and applications. We stress again that these mathematical 
issues do not render the MacGregor fc Charbonneaul 1 1997 ) analysis invalid or alter its conclusions; they are brought into 
stark light only when the new term arising from the tt x J effect is included. 



5 CONCLUSIONS 



In this paper, we test for the exist ence of interfacial dynamo in a protoneutron star. In Sec. 12.21 we apply a previous analysis 
I MacGregor fc Charbonneaulll997 ) pertaining to the situation where a and r\ are discontinuous, the discontinuities are located 
at the same place, and the disco ntinuities are located where the shear is also discontinuous. We confi rm that such a system 
never leads to a growing dynamo (jPeinzer fc StbJll97ll : lparkerll 19931 : IMacGregor fc Charbonneaulll997l ): when the a-effect and 
discontinuous velocity shear occur at the same location, the convective fluid motions tend to suppress the a-effect and dynamo 
amplification: the a-effect does not have enough time to generate the toroidal component of the field before being overcome 
by velocity shear ( Parker 19931 : MacGregor fc Charbonneau 1997 ). Note that there is currently no physical observation of a 
neutron star that can indicate whether the discontinuities in velocity, a, and r\ are coincident or not, unlike in the Sun, where 
helioseismology has revealed the existence, location, and the thickness of the tachocline (MacGregor & Charbonneau 1997). 

In Sec. 12.31 we set apart the velocity shear and a and r\ discontinuities by a distance d and apply the results to a 
protoneutron star. We find that there is no growth for k = 2kd > 1, basically because the discontinuities (where the field 
components are generated) get too far apart, and the waves damp before they can travel to the other layer; the velocity shear 
then cannot draw upon the toroidal component to regenerate the poloidal component, and the a-effect cannot draw upon the 
poloidal component to regenerate the toroidal component, before they diffuse away. As we increas e u, we find that the quickest 
growin g dynamo modes are those with k = 1, just like for the Solar dynamo case investigated by 
1997). At this value of k, we find that we require v w 195 to obtain growing modes, whereas 



MacGregor & Charbonneau 



MacGregor fc Charbonneau 



19971 ) found that they only need v > 20. Thus, we find that, for protoneutron star values of t) and a, the dynamo requires 



stronger shear and/or a-effect, compared to the Solar case, to be able to grow and be self-sustaining. This is because, in a 
protoneutron star, the upper region's magnetic diffusivity is smaller than the lower region's (the inverse of the Solar case), 
i.e. the waves dampen faster in the zone between the shear interface and the a interface, so the source terms at both ends of 
the dynamo process (proportional to 70, which controls the strength of the a-effect, and cjq, which controls the strength of 
the velocity shear) must be larger to compensate. 

As a specific application, we derive the condition for growing modes for a protoneutro n star (with typical values of 
77c = Wr]„f and R — 10 4 m) in Eq. (28). If we assume the values of a, AT?, and ljq given by Bonanno. Urpin. fc Belvedere! 
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(2005), we find that the threshold period for growing dynamo depends on R c , the location of the velocity shear: P < 49.8 ms 
if R c = 0.3R, P < 83.6 ms if R c = 0.6R. 

I n Sec. 3, we demonstrate that th e fix J dynamo effect, arising from anisotropic vorticity turbulence |Radlerlll969l ; lRobertsl 
19721 : iRoeachevskii fc Kleeorinll2003h . excites a dynamo in situations where it would otherwise be damped. For simplicity, 



we confine the effect to the boundary where the shear, a, and rj are also discontinuous. We find that the ft x J term in the 
induction equation can generate growing modes even when the discontinuities are coin cident, where one would previously 



expect zero growth |Deinzer fc Stixlll97ll ; |Parkerlll993l : iMacGregor fc Charbonneaull 19971 ). Designating the value of v\ [which 
characterizes the fi x J effect, defined in Eq. (43)] which gives growing modes for given n and v as zA. m m, we plot 

^l.min versus 

v for n — 10 and n = 3 in Fig. [4] We see in Fig. |4]that ^i, m i n approaches a constant value as v — » oo. Thus, we find that, for 
a particular n, there is a minimum value of v\ for which growth is guaranteed, regardless of v. We designate this minimum as 
v\ X and plot Vi^ versus n in Fig. [S] Applying these results to protoneutron stars, we derive the condition for growth due to 
the ft x J effect, with rj c = 10ry n /, in Eq. (44). 

In Sec. 4, we discuss so me physical and mathematical issu es regarding our calculations of the interfacial dynamo and 
those in the classic paper by IMacGregor fc Charbonneaul (|1997l ). Starting from the original induction equation, we find that 
the inclusion of a new term (arising from the ft x J effect) mean that there are some terms in v olving derivatives with respect 
to z which cannot be ignored in general but were neglected by MacGregor fc Charbonneaul (1993) an d ourselves earlier in 
this paper. We also find new problems with the boundary and jump conditions when this new effect is introduced: either the 
boundary and jump conditions are not enough to solve for all the unknowns (if we assume continuous and uniform fl x J 
effect) or they reduce the problem to trivial solutions (if we assume the CI x J effect is isolated at the interface, where the 
jump in fi is represented with a 5-function). We point out these fundamental issues without resolving them, a substantial task 
which lies beyond the scope of this paper. 

Throughout this paper, we consider only the dynamo operating near the interface between the convective and neutron 
finger regions of a protoneutron star. The reason for this is that the difference in diffusivity between the two regions is 
expected to make this dynamo more e fficient than the conventiona l a-fl dynamo. Without taking into account the different 
diffusivity in the neutron finge r region, iThompson fc Duncan! J 1993!) found that protoneutron star dynamos require P < 1 ms 
to operate efficiently; whereas Bonanno. Urpin. fc Belvedere! (|2005l ). who proposed a jump in diffusivity across the interface, 
fo und higher threshold periods for their dynamo, P < 1 s (i.e., the dynamo is easier to excite). This seems to be confirmed 
by IZhang. Liao. fc Schubertl (|2004l ) , who found that the generated fields are still c oncentrated in the tachocline r e gion, even 
if other layers away from it are included in the calculation. Furthermore, following iBonanno. Urpin. fc Belvedere ( 2005), we 



assum e that the dynamo operates in the period when the entire protoneutron star is unstable [< 40 s after its birth (|Naso et al 



2008)], i.e., there are no stable regions above the neutron finger region and below the convective region. Note, however, that 



the actual existence of the neutron finger instability r egion surrounding the convectively unstable region is still disputed 
( Buras et aDbood : Bessart et allbood : iNaso et al.ll2008r ). 
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